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^ . Abstract. This paper deals with the investigation of a closed form solution 

[^^ I of a generalized fractional react ion- diffusion equation. The solution of the 

'nI" I proposed problem is developed in a compact form in terms of the H-function 

Q ' by the application of direct and inverse Laplace and Fourier transforms. Frac- 

^ ■ tional order moments and the asymptotic expansion of the solution are also 

^ . obtained. 

^ ■ 

B '. 1 Introduction 



X 



It is a known fact that react ion- diffusion models play a very important role 
in pattern formation in biology, chemistry, and physics, see Wilhelmsson and 
Lazzaro (2001) and Frank (2005). These systems indicate that diffusion can 
produce the spontaneous formation of spatio-temporal patterns. For details, 
one can refer to the work of Nicolis and Prigogine (1977) and Haken (2004). 
A general model for reaction-diffusion systems is investigated by Henry and 
Wearne (2000, 2202) and Henry et al. (2005). 

The simplest reaction-diffusion models are of the form 

dN d'^N 



where d is the diffusion constant and F{N) is a nonhnear function rep- 
resenting reaction kinetics. It is interesting to observe that for F{N) = 
'~fV{l — N), eq.(l) reduces to Fisher-Kolmogorov equation and if we set 
F{N) = '~fN{l — N'^), it gives rise to the real Ginsburg-Landau equation. 
Del-Castillo-Negrete, Carreras and Lynch (2002) studied the front propaga- 
tion and segregation in a system of reaction-diffusion equations with cross- 
diffusion. Recently, del-Castillo-Negrete et al (2003) discussed the dynamics 
in reaction-diffusion systems with non-Gaussian diffusion caused by asym- 
metric Levy flights and solved the following model 

dN 

riD^N + F{N),N = N{x,t), (2) 



dt 
with F = 0. 

In this paper, we present a solution of a more general model of fractional 
reaction-diffusion systems (2) in which ^ has been replaced by the Riemann- 
Liouville fractional derivative o-Df , P > 0. The results derived are of a more 
general nature and than those investigated earlier by many authors, notably 
by Jespersen, Metzler, and Fogedby (1999) for anomalous diffusion and del- 
Castillo-Negrete et al. (2004) for the reaction-diffusion systems with Levy 
flights, and fractional diffusion equation by Kilbas et al. (2005). The solution 
has been developed in terms of the H-function in a compact and elegant form 
with the help of Laplace and Fourier transforms and their inverses. Most 
of the results obtained are in a form suitable for numerical computation. 
The present study is in continuation of our earlier works, Haubold (1998), 
Haubold and Mathai (2000), and Saxena, Mathai, and Haubold (2002, 2004a, 
2004b, 2004c, 2005). 

2 Mathematical Preliminaries 

A generalization of the Mittag-Leffler function (Mittag-Leffler, 1903, 1905), 

oo n 

^0 r(na + 1) 
was introduced by Wiman (1905) in the generalized form 

oo n 

^0 r(na + /3) 



The main results of these functions are available in the handbook of Erdelyi, 
Magnus, Oberhettinger, and Tricomi (1955, Section 18.1) and monographs 
by Dzherbashyan (1966, 1993). The H-function is defined by means of a 
Mellin-Barnes type integral in the following manner (Mathai and Saxena, 
1978), 
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and an empty product is always interpreted as unity; tti, n,p,q G Nq with 
Q <n <p,l <m < q, Aj,Bj E -R+, Oj, hj G i? or C{i = 1, . . . ,p;j = 1, . . . , g) 
such that 

Ai{bj + k) ^ Bj{ai - I - l),k,l e No;i = 1, . . . ,n;j = 1, . . . ,771, (7) 

where we employ the usual notations: A^o = (0, 1, 2, . . . , ); i? = (— oo, oo), i?+ = 
(0, oo) and C being the complex number field. The contour fi is either 
L_oo) -^+00 or Li^oo- The following are the definitions of these contours. 

(i) Q = L_oo is a left loop situated in a horizontal strip starting at the 
point —00 + iipi and terminating at the point —00 + iip2 with —00 <(/?!< 

(P2 < +00; 

(ii) Q = L+00 is a right loop situated in a horizontal strip starting at the 
point +00 + iipi and terminating at the point +00 + iip2 with —00 < ipi < 

(P2 < +00. 

(iii) Q = Li^oD is a contour starting at the point 7 — ioo and terminating 
at the point 7 + ioo, where 7 G -R = (—00, +00). 

A detailed and comprehensive account of the H-function is available from 
the monograph by Mathai and Saxena (1978), Prudnikov et al. (1989), and 
Kilbas and Saigo (2004). The relation connecting the p\l'g(z) and the H- 
function is given for the first time in the monograph by Mathai and Saxena 



1978, p.ll, Eq.1.7.8) as 

,Tr r(ai>-4i),...,(ap,Ap)|i 
P^qi(bi,Bi),...,(B„B,)\^\ 






(l-ai,yli),...,(l-ap,Ap) 1 

(0,l),(l-6i,Bi),...,(l-b„fe5)J' 



where p\E' 5(2;) is Wright's generalized hypergeometric function, (Wright, (1935, 
1940); also see Erdelyi, Magnus, Oberhettinger, and Tricomi (1953, Section 
4.1), defined by means of the series representation in the form 



P^g(^) 



-P ^ikbq.Bq) \^\ 



E 



^,[n]^,T{h,+Bf)]{r)V 



(9) 



where z & C, ai, bj G C, Ai, Bj G -R+; 

X]J=i &j — Yl'j=i^j > —^',C being the complex number field. This function 

includes many special functions besides the Mittag-Leffler function defined 

by the equations (3) and (4). It is interesting to observe that for Ai = Bj = 

1 \/ i and j, (9) reduces to a generalized hypergeometric function pFq{z). 

Thus 






T^un 
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(10) 



where aj 7^ — z/ {j = 1, . . . ,p and z/ = 0, 1, 2, . . .);p < g (or p = g, |z| < 1). 
Prior to (9), Wright (1933) introduced a special case of (9) in the form 



<^{a,b;z) = o^i[ 



{b,a) I 



E 



g r(ar + b) (r)!' 



which widely occurs in problems of fractional diffusion. It has been shown 
by Saxena, Mathai, and Haubold (2004b) that 



E. 



a,l3\ 
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[(/3,a)l 



mii-^ 



(0,1) 1 

(0,l)(l-/3,a)J- 



If we further take /3 = 1 in (11) and (12), we find that 

E^,,iz) = E^{z)= i^i[[S)k] 



Hl^^i-z 



(1,1) 1 

(0,1) 
(0,1), (0,a)- 



(11) 

(12) 

(13) 
(14) 



where Re{a) > 0, a G C. The definitions of the well-known Laplace and 
Fourier transforms of a function N{x, t) and their inverses are described 
below. 

The Laplace transform of a function N{x, t) with respect to t is defined 

by 



L{N{x,t)} 



e'''*N{x,t)dt, t>0,xeR, 



(15) 



where Re{s) > 0, and its inverse transform with respect to s is given by 

\ ff+iOD 



L-^ {N{x, t)} = -^ I e''N{x, s)ds, 

2m J-f-ioo 



(16) 



7 being a fixed real number. 

The Fourier transform of a function N{x, t) with respect to x is defined 

by 



F{N{x,t)} 



e'^''N{k,t)dk. 



The inverse Fourier transform with respect to k is given by the eq. 

1 /-oo 



F-'{N{x,t)} = — I e 

ZTT J -oo 



-ikx 



N{k,t)dk. 



(17) 



(18) 



The space of functions for which the transforms defined by (15) and (17) 
exist is denoted by LF = L{R+) x F{R). 

In view of the results of Saxena, Mathai, and Haubold (2004a, p. 49), also 
see Prudnikov, Brychkov, and Marichev (1989, p. 355, Eq.2.25.3.2), the cosine 
transform of the H-function is given by 



/o" tP-^cos{kt)H!;^g''[ar \^^^^^^\bg, Bg)]dt 

"/C^ (l-fe,,iJ,),(l±£,|) 

(p,A.),(l-ap,Ap),(i±£,f) 



(19) 






where Re[p + fxY<f<n.{^)] > 1, \arg a\ < ivr^; 9 = E"=i A, - E^=„+i Aj + 

The Riemann-Liouville fractional integral of order u is defined by (Miller 
and Ross, 1993, p. 45) 



oDrN{x,t) 



T(u) Jo 



{t-uy-^N{x,u)du, 



(20) 



where Re{u) > 0. 

Following Samko, Kilbas, and Marichev (1990, p. 37), we define the 
fractional derivative of order a > in the form 

^n^T, N 1 (i" /■* N(x,u)du „ , , , _ ,„_, 

°^°^'^-'> = f(^r37^5Fi ft-V-... ^ t>0.(n = |a.] + l. (21) 

where [a] means the integral part of the number a. 

A comprehensive account of the various applications of fractional calculus 
in physics is available from the monograph of Frank (2005). 

From Erdelyi et al. (1954a and 1954b, p. 182), we have 

L{oD;''N{x,t)} = s''F{x,s), (22) 

where F{x, s) is the Laplace transform with respect to t of N{x, t), Re{s) > 
and Re{u) > 0. 

The Laplace transform of the fractional derivative, defined by (21), is 
given by Oldham and Spanier (1974, p. 134, Eq.(8.1.3)) 

n 

L{oD^N{x,t)} = s''N{x,s)-J2s'''^ oD^-''N{x,t)\t=o,n-l < a < n. (23) 

r=l 

In certain boundary value problems the following fractional derivative of 
order a is introduced by Caputo (1969) in the form 

oDtfix,t) = — ^ r /^""H^'^)^^ (24) 

* •^ ^ ' ^ T{m - a) Jo {t - r)°+i-'" ^ ^ 

m — 1 < a < m, Re{a) > 0, m ^ N, 
. ?!|(^.U„ = », (25) 

where -^f is the m^^ derivative of order m of the function f{x,t) with 
respect to t. The Laplace transform of this derivative is given by Podlubny 
(1999) in the form 

m— 1 

L{QD^f{x,t);s} = s''F{x,s)- J] s''~'''^f'^''\x,0+),m-l< a < m. (26) 

r=0 

The above formula is useful in deriving the solution of differential and integral 
equations of fractional order governing certain physical problems of reaction 
and diffusion . 
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We also need the Weyl fractional derivative, defined by 

D^fix t) - ^ ^ r ^^^'^^^^ (27) 

where x E R, fi > 0,n = [fi] + l,[fi] being the integer part of /i > (Samko 

et al, 1990, Section 24.2). 

Its Fourier transform is given by (Metzler and Klafter, 2000, p. 59, A. 11) 

F {_oo D^fix, t)} = {ikr^ik, t), /i > 0, (28) 

where \E'(/c, t) is the Fourier transform of /(x, t) with respect to the variable x 
of f{x,t). Following the convention initiated by Compte (1996), we suppress 
the imaginary unit in Fourier space by adopting the slightly modified form 
of above result in our investigations(Metzler and Klafter, 2000, p. 59, A. 12) 

F{_^D^J{x,t)} = ~\kr^{k,t) (29) 

instead of (28). Finally we also need the following property of the H-function 
(Mathai and Saxena 1978) 



where 6 > 0. 



TTm,nt 5 '"- ''"^■' 



(ap,Ap)i 



t^pT[^ (w;a (30) 



3 Fractional React ion- Diffusion Equation 

In this section, we will investigate the solution of the generalized reaction- 
diffusion equation (31). The result is given in the form of the following 
Theorem. Consider the generalized fractional-reaction diffusion model 

oD^N{x, t) = r]_^D^N{x, t) + v?(x, t); ?] > 0, t > 0, a; G i?, 1 < /9 < 2, (31) 

< q; < 1 with the initial conditions 

^Dt-^N{x,Q) = f{x),oDr^N{x,0) = g{x),xe R,lim^^±^N{x,t) = 0, 

(32) 
where o-D" ^iV(x, 0) means the Riemann-Liouville fractional derivative of 
order a — 1 with respect to t evaluated at t = 0. Similarly, qD'^^'^N{x, 0) 
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means the Riemann-Liouville fractional derivative of order a — 1 with respect 
to t evaluated at t = 0. The quantity 77 is a diffusion constant and v^(a;, t) 
is a nonlinear function belonging to the area of reaction-diffusion. Then for 
the solution of (31), subject to the initial conditions (32), there holds the eq. 



-1 



/;/3— i /-OO 

ZTT J-00 

+ -TT- / ~g{k)Efi^fi^x{-r]\k\°'t^)exv{-ikx)dk 
Zn J-00 

+ 7^ / i^~' / V^{k,t-i)Ep^p{-7^\k\^e)exp{-ikx)dkdi, 
Zn JO J-00 

where ~ indicates the Fourier transform with respect to the space variable 

X. 

Proof. If we apply the Laplace transform with respect to the time variable 
t and use eq. (23), eq. (31) becomes 

s^N*{x, s) - fix) - sg{x) = r]^^D^N*{x, s) + ^*{x, s). 

As is customary, it is convenient to employ the symbol * to indicate the 
Laplace transform with respect to the variable t. 

Now we apply the Fourier transform with respect to the space variable x 
to the above equation and use the initial conditions and the result (29), then 
the above equation transforms into the form 

Jy [k, S) = —. —— + — r —— + — r —— , (34) 

where, according to convention followed, ~ indicates the Fourier transform 
with respect to x. On taking the inverse Laplace transform of (34) and using 
the result 

^'' l^r-^'4 = i""^^a,a-/3+i(-at"), (35) 

where Re{^s) > 0, Re{a — jS) > —1, it is seen that 

N{k,t) = f{k)t''-'Ep^p{-r]\krt'') + ~g{k)t^~'Ep,p^ii-v\krt'') 

+ U{k,t-i)e''Ep,p{-^\k\^e)di. (36) 

JO 



The required solution (33) now readily follows by taking the inverse Fourier 
transform of (36). Thus, we obtain 

N{x,t) = —— / f{k)Ep,p{-r]\k\''t^)exp{-ikx)dk 



2n 


J — CO 


^f3-2 


(•OO 


2n 


i-oo 


1 / 


rt 



g{k)Ep^p_i{-r]\k\''t'')exp{-ikx)dk 

1^ ff'' I V*ik,t-OEpA-v\kn^)expi-tkx)dkd^. 
zvr Jo J-oo 

This completes the proof of the theorem. 

Note 1. It may be noted here that by virtue of the identity (12), the solu- 
tion (33) can be expressed in terms of the H-function as can be seen from 
the solutions given in the special cases of the theorem in the next section. 
Further we observe that (33) is not an explicit solution and special cases are 
of interest. 



4 Special Cases 

When g{x) = 0, then applying the convolution theorem of the Fourier trans- 
form to the solution (33), the theorem yields 
Corollary 1.1. The solution of fractional reaction-diffusion equation 

oD^N{x, t) = ri_^D^Nix, t) + ^{x, t), t > 0, r^ > 0, (37) 

subject to the initial conditions, 

^Dr^N{x,t)\t=, = f{x),oDr'N{xM=o = (38) 

for X e R, limx^±aoN{x, t) = 0,1 < /3 <2, 

^ a ^ 1, where 77 is a diffusion, constant and ^p{x, t) is a nonlinear function 
belonging to the area of reaction diffusion is given by 



N{x,t) = / Gi{x-T,t)f{T)dT 
J —00 

+ f\t - 0^-' r G2{x -r,t- O^r, Odrdi, (39) 

JQ Jo 



where 



Gi{x,t) 



t" 



-1 



27r j-oo 



exp{-ikx)Ep^p{-r]\k\'^t^)dk 



t^'^ 2 1 

a 



/ cos{kx)Hl'l[kri^H' 
Jo 



l/a.p/a 



\X\ 



X\ 3,3L^i/„^/3/a 



(l,l/a),(/3,/3/a),{l,l/2) 
(l,l),(l,l/a),{l,l/2) 



(0,1/a) -1 ,, 

(0,l/a),(l-/3,/3/a)J«'« 

, i?e(a) > 0, (40) 



G2{x,t) 



1 /""^ 
2n J-oo 



exp{—ikx)Ep^p{—rj\k\°'t^)dk 



— / cos(kx)Hl'l[k7]^/''t^/'' 
ira Jo 



(0,1/a) 1 .. 

(0,l/a),{l-/3,/3/a)J"'« 



r//- 



2,1 r 



\X\ 



a 



X\ 3.3l^l/a^/3/a 



(l,l/a),(/3,/3/a),{l,l/2)n 
(l,l),(l,l/a),{l,l/2) J' 



Re{a) > 0. (41) 



If we set /(x) = S{x), where 6{x) is the Dirac delta function, then we arrive 

at 

Corollary 1.2. Consider the following reaction-diffusion model 



d/^Njx, t) 



r] _^D'^N{x,t), r7>0, x e R, 



(42) 



with the initial condition N{x,t = 0) = 6{x),limj.^±ooN{x,t) = 0,0 < /3 < 
1, where rj is a diffusion constant and 6{x) is the Dirac delta function. Then 
the solution of (42) under the given initial conditions is given by 



N{x,t) 



tC 



-1 



a\x\ 



:H. 



2,lr 
3,3 [ 



\X\ 



[r^t^y/c 



{l,l/a),(/3,/3/a),{l,l/2). 
{l,l),(l,l/a),{l,l/2) h 



(43) 



where Re{a) > 0, i?e(/3) > 2. 

When f3 = 1/2, the above corollary reduces to the following result. 
Consider the following reaction-diffusion model 



d^/^N{x,t) 



T] ^^D'^N{x,t), 77>0, xeR, 



(44) 



with the initial condition N{x,t = 0) = 6{x),lim.x^±ooN{x,t) = 0, where r) 
is a diffusion constant and 6{x) is the Dirac delta function. Then the solution 
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of (42) under the given initial conditions is given by 



Nix,t) = —-^H, 



1 --2,lr fI 



a|x|tl/2 3,3L^^^i/2)l/c 



(1, l/a),{l/2,l/2a), (1,1/2)- 
(l,l),{l,l/a),{l,l/2) 



where /2e(a) > 0, Re{(3) > 0. 

Remark 1. The solution of eq. (42), as given by Kilbas et al (2005) is 
in terms of the inverse Laplace and inverse Fourier transforms of certain 
functions whereas our solution of the same equation is obtained in an explicit 
closed form in terms of the H-function. 

An interesting case is, when /5 — > 1, then in view of the cancellation law 
for the H-function (Mathai and Saxena, 1978), (43) provides the following 
result given by Jespersen et al. (1999) and recently by del-Castillo- Negrete 
et al. (2003) in an entirely different form. 

For the solution of fractional-reaction-diffusion equation 

^^N{x,t)=r]_^D:N{x,t), (45) 

with initial condition 

N{x,t = 0) = S{x), limx^±ooN{x,t) = 0, 
there holds the relation 



1 „iir la;| 



a\x\ 2,2L^i/q^i/„ 



(1,1/a), (1,1/2)- 



N{x, t) = -^Mi [-rBrr. r^V^^ir ] > (46) 



where Re{a) > 0. In passing, it may be noted that (46) is a closed form 
representation of a Levy stable law, see Metzler and Klafter (2000). It is 
interesting to note that as a ^ 2, the classical Gaussian solution is recovered 



as 



N(r f) - J-H^'^\^^ (i,i/2),(i,i/2) 
^'J[x,t) - 2\x\'^'^\r]ty-^ (i-i)-(i-i/2) 



2\x\ '''\r]ty/^ (I'l) 

(47rr7t)-^/2exp(-M). (47) 
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It is useful to study the solution (43) due to its occurrence in certain 
fractional diffusion models, we will find the fractional order moments of (43) 
in the next section. 

Remark 2. Applying Fourier transform with respect to x in (42), it is found 
that 

—^{k,t) = -7^\kr^{k,t), 0<f3<l, (48) 

which is the generalized Fourier transformed diffusion equation, since for 
a = 2 and for /? ^ 1 , it reduces to the Fourier transformed diffusion equation 

being a relaxation (diffusion) equation, for a fixed wave number k (Metzler 
and Klafter, 2000). Here \E'(/c, t) is the Fourier transform of N{x,t) with 
respect to x. 

Remark 3. It is interesting to observe that the method employed for de- 
riving the solution of the equations (31) and (32) in the space = LF = 
L{R^) X F{R) can also be applied in the space LF' = L'{R^) x F', where 
F' = F'{r) is the space of Fourier transform of generalized functions. As 
an illustration, we can choose F' = S' or F' = D'. The Fourier transforms 
in S' and D' are introduced by Gelfand and Shilov (1964). S' is the dual 
of the space S , which is the space of all infinitely differentiable functions 
which, together with their derivatives approach zero more rapidly than any 
power of l/\x\ as \x\ -^ oo (Gelfand and Shilov, 1964, p. 16). D' is the dual 
of the space D which consists of all smooth functions with compact supports 
(Brychkov and Prudnikov, 1989 p. 3). For further details , the reader is 
referred to the monographs by Gelfand and Shilov (1964) and Brychkov and 
Prudnikov (1989), if we replace the Laplace and Fourier transforms in (15) 
and (17) by the corresponding Laplace and Fourier transform of generalized 
functions. 

5 Fractional Order Moments 

In this section, we will calculate the fractional order moments, defined by 

< \x\^ >= r \x\^N{x,t)dx. (49) 
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Using the definition of tlie Mellin transform 

M{f{ty,s)}= / f-'f{t)dt, (50) 

JO 

we find from (49) that 

/oo 
\x\^N{x,t)dx (51) 

-oo 

< \x\^(t) >- ?^ ra;'5-lH-2,lr^^_ (l,l/a),(l,/3/a),(l,l/2). , .52) 

< Fl W >- ^ 7o ^ ^3,3[^i/„^/3/c, (l,l),(l,l/a),(l,l/2) J"^' l^^J 

Applying the Melhn transform formula for the H-function, namely 

j^x^-'H^^;^-[ax \%lf:^]dx = a-^e(-p), (53) 

where 



|) < i?e(p) <r<7 <„ i?e(l-^"- ' ' ^ 



r<7<™^e(-^) < Re{p) <r<7<„ i?e(-^), |ar^ a\ < -n9, 9 > 0, 



is defined in (19) and 0(— p) in the definition of the H-function (5), we see 
that 



2 A/a /3{(5/q)+i_(i//3)} r( ^)r(i + ^)r(i + ^ , 
r(-|)r(/3 + f)r(i + 



< |a;|^(t) >= ±r^s/a^m/a)+i-(m} ' ^ J' ^' ^'' ^' + J (54) 



whenever the gammas exist, Re{6) > —1 and Re{6 + a) > 0. 
Two interesting special cases of (54) are worth mentioning. 

(i) As 6 ^ 0,then using the result ~ 2; for 2; << 1, (55) 

1 [z) 

we find that 

^-:o < \x\\t) >= ^. (56) 

(ii) When a = 2,S = 2, the linear time dependence, 

of the mean squared displacement is recovered. 
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5_2,.-.2<kWr>=-F7TT:5r' (57) 



6 Behavior of the Solution in Equation (43) 

Eq. (43) can be expressed in terms of the Mellin-Barnes type integral (Erdelyi 
et al, 1953, chapter 1) as 



N{x,t) 



t^'^ 1 



a\ 



r(^)r(i-s)r(i 



\x\ 



x\ 27r Jl r(/5 - f )r(l - f )r(s/2) ^r/V"t/3/^ 



'ds. 



(5J 



Let us assume that the poles of the integrand of (58) are simple. Now 
evaluating the sum of residues in ascending powers of |x| by calculating the 
residues at the poles of r(l — s) at the points s = 1 + z/ (z/ = 0, 1, 2, . . .), and 
r(l — -) at the points s = 1 + z/ (z/ = 0, 1, 2, . . .), it is found that the series 
expansion of the general solution (43) is given by 



N(x,t) 



t'-' ^ {-ly r[i±^]r(i - (1 + u)/a) 



arj^lnt^l' 



u=0 



\x\ 



^l/a-j-P/c 



+ 



\X 



a— I oo 



•qt 



■E 



(-irr(i-a(i + z/)) 
r(-/5z/)r(i-f(i + z/))r(^^) [{ritPY/-_ 



(59) 



where < Re{v) < a, { ^i/L'f'i/c } < 1- 
From (59), we infer that 



N{x, t) ~ At^-^/°-i + 5|x|""\ as X -^ 0, 



(60) 



where A and B are numerical constants . Further from the series expansion 
(59), it can be seen that the initial behavior is given by 



N{x,t) 



r(i + i)r(i-i) 



for 1 < a < 2, 



(61) 



with {j:;^} « 1. 

Next, if we calculate the residues at the poles of r{s/a) at the points 

s = —av (z/ = 0, 1, 2, . . .), it gives 



N{x,t) 



t" 



— 1 oo 



E 



Til + av) 



\x\ — r(/^ + /3z/)r(i + ^)r( 



u=0 



2 > 



r]t^ 



\x\ 



(62) 
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for |a;| > 1. From (62), it can be readily seen that 

N{x,t) 



\x\ 



for large |x|. In conclusion, it is observed that the solution given by (43) does 
not admit a probabilistic interpretation in contrast with fractional reaction- 
diffusion based on Caputo derivative derived by the authors. However, when 
13 ^ 1 then it has a probabilistic interpretation, as shown in corollary 1.2. 
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